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Quantum manipulation protocols for quantum sensors and quantum computation often require 
many single qubit rotations. However, the impact of phase noise in the field that performs the 
qubit rotations is often neglected or treated only for special cases. We present a general framework 
for calculating the impact of phase noise on the state of a qubit, as described by its equivalent 
Bloch vector. The analysis applies to any Bloch vector orientation, and any rotation axis azimuthal 
angle for both a single pulse, and pulse sequences. Experimental examples are presented for several 
special cases. We apply the analysis to commonly used composite 7r-pulse sequences: CORPSE, 
SCROFULOUS, and BBl, used to suppress static amplitude and detuning errors, and also to spin 
echo sequences. We expect the formalism presented will help guide the development and evaluation 
of future quantum manipulation protocols. 
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I. INTRODUCTION 

Atomic quantum sensors and tests of fundamental 
physics commonly rely on the ability to rotate a Bloch 
vector representing a spin-i system or qubit. Besides 
quantum state manipulation, rotations can also be used 
to undo inhomogeneous errors or to reduce other sources 
of noise, as is done with spin echo pulses or dynamical 
decoupling [IHZ|- Precise rotations, required for manip- 
ulating collective spin-squeezed states [8Hl6]. might also 
enable dynamical spin-squeezing to the two-axis squeez- 
ing limit |17i using already realized one-axis twisting in 
cold atoms systems [TMT?1 Hi] . Given that actual rota- 
tions are imperfect, it is an open question whether such 
rotation protocols can be realized without adding large 
amounts of additional noise, thus destroying the squeez- 
ing. 

Most rotation protocols assume that the phase of 
the field that rotates the qubit is perfectly stable, and 
that imperfections arise only due to slowly varying am- 
plitude errors or detuning errors. Composite rota- 
tions sequences [T51 - E^ and generalizations to shaped 
pulses [131 m], including optimal control theory 25 27J, 
can be used to reduce these errors to essentially arbitrary 
order. 

In reality, the phase of the qubit-field coupling is never 
perfectly stable, largely due to phase noise in the local 
oscillator (LO) used to generate the field. The LO is 
typically a radio or microwave oscillator in nuclear spin, 
superconducting Josephson junction, quantum dot, neu- 
tral atom/ion Zeeman and hyperfine qubit systems |28j . 
In the case of highly-forbidden optical transitions [53], 
the LO is an ultra-stable laser. Further, qubit transition 
frequency fluctuations can be straightforwardly mapped 
onto an equivalent phase noise of the LO. Such a fluctua- 
tion might arise due to noise in the DC bias current of a 
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superconducting Josephson junction qubit [30j . or differ- 
ential light shifts for atomic qubits in an optical trap [29j . 
It is critical for future work beyond proof-of-principle ex- 
periments to develop general tools for analyzing the im- 
pact of phase noise on a rotation, both in terms of over- 
all fidelity for quantum gates, and quadrature specific 
noise for manipulating states with anisotropic sensitiv- 
ity to noise, such as spin-squeezed states [SHU], reduced 
spin- noise states [31], and Dicke states [32] . 

The effect of phase noise on atomic response has been 
studied in various contexts, for example on atomic exci- 
tation probability [33 - 36] under continuous drive. In the 
atomic sensor community, the impact of phase noise is 
evaluated for a single quadrature for specific Ramsey se- 
quences [57H55] . In contrast, we present a general frame- 
work that may be applied to an arbitrary LO phase noise 
spectrum for continuous resonant drive without making 
assumptions about the orientation of the Bloch vector. 
The framework can be extended to arbitrary resonant 
pulse sequences assuming white LO phase noise. Fur- 
thermore, we fully specify all second order noise moments 
of the Bloch vector including covariances and variances, 
important for predicting the fidelity of single qubit gates 
and manipulations of spin-squeezed or Dicke states. The 
methodology and tools presented in this paper can help 
guide the development and evaluation of future quantum 
control and measurement protocols. 

The organization of the paper is as follows. In Sec. |ITJ 
we describe the qubit-field interaction as a rotation of 
a Bloch vector. We show that the net effect of a co- 
herently phase modulated rotation can be reduced to a 
small rigid rotation of the Bloch sphere. A description 
of the experimental system used to demonstrate the the- 
ory is provided in Sec. |I11| Experimental examples of the 
response of the Bloch vector driven by phase modulated 
rotations are then presented. In Sec. |IV[ we extend the 
analysis in Sec. [IT] assuming linear response to relate the 
single sideband (SSB) LO phase noise to Bloch vector 
noise projections through a covariance transfer matrix 
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for a single rotation. Experimental realizations for a few 
special cases are also presented. The Bloch vector noise 
projection variances and covariances are captured in the 
covariance noise matrix. In Sec. |V[ we generalize the 
single-rotation covariance noise matrix to that for mul- 
tiple rotations/pulse sequences and apply the tools to 
commonly used composite 7r-pulse sequences: CORPSE, 
SCROFULOUS, and BBl [18..22j, and also to spin echo 
pulse sequences. A simple formula for the average infi- 
delity of any pulse sequence, appropriate in the context 
of quantum computing, is also presented. Finally, we give 
a summary of the results presented in this paper, and an 
outlook of future work in Sec. IVII 



II. DEFLECTION OF BLOCH VECTOR DUE 
TO COHERENTLY PHASE MODULATED 
ROTATION 

A. System Hamiltonian 




The system considered in this paper consists of a qubit 
with transition frequency fa coupled to a resonant clas- 
sical electromagnetic field with strength characterized by 
the Rabi frequency fn proportional to the field ampliu- 
tude. The system's Hamiltonian in the lab frame is 

-fflab = hfaCT:, + hfn COS (27r/at + , (1) 

where h is the Planck constant, (f){t) is the LO phase as 
a function of time and ,<Jy, are the Pauli matri- 
ces. The first term corresponds to the Hamiltonian of the 
qubit in the absence of the field, and the second term cor- 
responds to the qubit-field interaction that drives Rabi 
flopping between the qubit states |t) and \\). By going 
into a rotating frame at the qubit transition frequency, 
and making the rotating wave approximation, the dy- 
namics is described by the Hamiltonian 

H,ot = ^ (cos(0(t))(T, -f sin(0(i)) ay) . (2) 

B. LO Phase Noise 

The LO phase (j){t) is random as a function of time, 
and its statistical properties can be captured by the au- 
tocorrelation function {(f>(t)(l>(t + T))t in the time domain, 
or equivalently the power spectral density of phase fluc- 
tuations S^{fm) in the frequency domain, where /,„ is 
the frequency offset from the LO carrier frequency Jlo- 
The power spectral density of phase fluctuations S^{fm) 
is related to the autocorrelation function through the 
Wiener-Khinchin theorem (see Appendix [A|) . The SSB 
phase noise £{fm) = S^{fm)/'2 of a LO is usually spec- 
ified in manufacturer datasheets instead of Scf){fm)- To 
facilitate application of our results in an experimental 
context, the SSB phase noise C{fm) is used throughout 
this paper. 



FIG. 1. (color online). Bloch Sphere Picture, (a) Geo- 
metrical representation of the Bloch vector and a coherently 
phase modulated rotation on the Bloch sphere. Initial Bloch 
vector Ji is specified by polar angle 8i and azimuthal angle 
In the absence of phase modulation, the rotation axis is at 
a constant angle (j>R from the x-axis. With phase modulation, 
the rotation axis oscillates in the x — y plane with amplitude P 
around its average position, (b) Elfect of a coherently phase 
modulated rotation. An unmodulated rotation R[(j>R,ijj) ro- 
tates Ji to J^. With phase modulation of the rotation axis, 
the final Bloch vector J/ is deflected slightly by the vector j/ 
that can be regarded as arising from a small rotation r. 

The phase noise spectrum of a LO is typically 
parametrized as a sum of 1// type phase noise as fol- 
lows 

oo „ 

where the coefficient Ck characterizes the strength of the 
noise spectrum with dependence 1//^- An important ex- 
ample is white phase noise where there is no frequency de- 
pendence. In this case, the only non-zero coefficient is for 
k — Q corresponding to C{fm) = C,o, and the autocorre- 
lation function is given by {(j){t)(j){t + T))t = CoS{t) where 
(5(t) is the delta function. Other frequently encountered 
phase noise spectra include flicker noise C{fm) — l//m 
and phase diffusion noise C{fm) = ^/ Im- 

C. Bloch Sphere Picture 

The evolution of the qubit's state in the rotating frame 
is conveniently visualized using the Bloch sphere picture. 
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A qubit's state can be mapped onto a Bloch column vec- 
tor J with components Jk = (cfc), where k = {x,y,z\. 
The tip of the Bloch vector resides on a Bloch sphere with 
radius equal to 1 (see Fig. [l^). The LO phase sets 
the instantaneous rotation axis around which the Bloch 
vector rigidly rotates. For fixed LO phase </)(t) — (pn, an 
initial Bloch vector Jj is mapped to an ideal final Bloch 



vector J J as 



R{(l)R,'iJj)3i 



(4) 



where R{4>r, ip) performs a counterclockwise rotation 
through an angle ip > about an axis in the x — y plane 
with azimuthal angle measured relative to x. The 
rotation angle = 2Tr fat is set by the Rabi frequency 
ffj and the amount of time t the field is applied. 



where labels on r have been suppressed. The relation- 
ships between the vectors Ji, J^, J/,j/, and the rotations 
R{(f)ii,^lj), r are depicted in Fig.[T]3. 

Because the dynamics with modulation can be de- 
scribed in terms of infinitesimal rigid rotations, the whole 
Bloch sphere is rigidly rotated by this additional phase 
modulation contribution, and r does not depend on J^. 
This is confirmed by solving for r analytically. Simple 
rotations about z through angle cjjj^ relate the small ro- 
tation r evaluated at 0^ = to the small rotation eval- 
uated at arbitrary 0/? through 

^ r{<l>R) = i?.(0fl)KO)i?._(-0fl) . (7) 
Therefore, it is sufficient to determine r for the special 
case (jjR — and 3i = x. Hereafter, quantities with a 
tilde (^) overhead apply to this special case only. We 
find 



D. Small Rotation describing Phase Modulated 
Rotation 

Our goal is to generate a general covariance trans- 
fer matrix that maps LO phase noise onto noise pro- 
jections of the final Bloch vector. To do so, we con- 
sider here the response of the Bloch vector to a co- 
herent phase modulation of the LO phase in the fre- 
quency domain, and generalize to a phase noise process 
in Sec. |IV[ A sinusoidal modulation of the LO phase 
(f){t) = + /? sin(27r/„ii -I- am) represents an oscilla- 
tion of the instantaneous rotation axis in the x — y plane 
with amplitude /3 about its average azimuthal angle (j)^ 
as illustrated in Fig. [T^. The modulation frequency and 
phase is fm and a^n respectively. 

A phase modulated rotation causes a defiection of the 
final vector J f from its ideal final orientation without 
modulation by 



(5) 



To first order in the phase modulation amplitude /3 <C 1, 
the deflection is perpendicular to Jy. The total effect of 
the modulation can be described by an additional small 
rotation r(0_R, ip, (3, /„, a™) applied after the ideal rota- 
tion such that 



"(0) = Ry{~J,)R,{jy) 




(8) 



where j/ = iO,jy,jz), and only first order in the small 
quantities jy and jz are retained. 

E. Solution for Deflection Vector 

We now solve for the jy and jz that determine the 
small rotation matrix r defined by Eq. [7] and [8] Writing 
the Heisenberg equations of motion for the Bloch vector 
components yields to first order in the small modulation 
amplitude 



dj± 

dip 



+ ij±_ = —/3 sin {xip + am) , 



(9) 



rJ}, 



(6) 



where j± = jz +'ijy, and x — fm/ fn- The rotating wave 
approximation fm <SC fa has been made. Note that this 
is the equation of motion for the coupled position and 
momentum of an undamped simple harmonic oscillator 
with natural resonance frequency fn driven with an ex- 
ternally applied force at frequency fm- The harmonic 
oscillator's displacement and velocity map onto jy and 
—jz respectively. Solving for jy and jz using the initial 
condition j_\_{ip = 0) = 0, we obtain 



Jy 

jz 



—x cos Um sin "0 — sin am cos ip + sin [xip -(- a™) 
x cos am cos Ip — sin am sin ip — x cos {xip + am ) 



(10) 
(11) 



The above solutions can be understood as a superpo- driven harmonic oscillator. The terms proportional to 
sition of the "transient" and steady state response of a 
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s'mip and cosip in Eq. 10 and 11 correspond to the re- 



(a) 



sponse of the harmonic osciUator at its natural frequency 
/fl. This response is called the transient response in 
damped harmonic oscillator systems. Because there is no 
damping in this oscillator, the "transient" response does 
not decay away. The terms proportional to sin {x^ + am) 
and cos (xip + am) correspond to the steady state re- 
sponse of the oscillator at the drive frequency fm- At 
X = 1, corresponding to the case of driving on resonance, 
the solutions take on the following limits 



\im jy 



1 



- cos {tjj + am) — cos a„i sin ip) , (12) 



lim = - - {ip sin {ip + am) + sin am sin tJj) . (13) 

The amplitude of the response grows roughly linearly 
with tfj for large ip ^ 1 a.s the drive is phase coherently 
adding momentum to the oscillator. 



III. EXPERIMENTAL DEMONSTRATION 

In order to connect the theory described in Sec.|IT]to an 
actual physical system, we experimentally demonstrate 
the linear response of the Bloch vector to a coherently 
phase modulated rotation for a few special cases. 

A. Physical Implementation 

The experimental system used for these studies was 
used to generate conditionally spin-squeezed states and is 
described in Ref. ^13] . We use an ensemble of = 7 x 10^ 
^''Rb atoms laser cooled and trapped in an optical lattice 
at 823 nm (Fig. The clock states |F = 2,mj^ = 0) = 
It) and \F — l,nip — 0) = ||) constitute a pseudo spin-i 
system/qubit with transition frequency fa = 6.834 GHz. 
All atoms are initially optically pumped into \\.). Follow- 
ing that, a microwave 7r/2-pulse rotates the Bloch vector 
up to the equator, initializing the system for the experi- 
ments. To a very good approximation, the effects of mi- 
crowave amplitude and phase inhomogeneity across the 
atomic ensemble may be neglected in our experiments. 
The intrinsic phase noise of the microwave LO, Agilent 
E8257D, is sufficiently low that we could use a small mod- 
ulation amplitude (3 in order to remain in the linear re- 
sponse regime and yet not be affected by phase noise of 
the LO source. Details on how phase noise from the mi- 
crowave LO affects measurement signal to noise can be 
found in [iO] . 

An experiment typically consists of a test rotation Rt 
using phase modulated resonant microwaves coupling the 
two-level system with Rabi frequency fa = 40.4 kHz. Af- 
ter the rotation is completed, the Bloch vector projection 
Jz = (A^t ~ ^i)/{^t + ^i) is obtained using a cavity- 
aided nondemoliton measurement of the state popula- 

tions iV^^^ = E-Li(l ± (^^))/2, where al is the Pauli 
spin operator corresponding to spin i. The deflection jz 
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FIG. 2. (color online). Experiment Schematic, (a) The 

clock states |t) and |^) of *^Rb form a pseudo spin-| system. 
Coupling of the atoms to a cavity mode resonant with the 
t) to je) transition creates a collective vacuum Rabi split- 
ting Qf which is probed to deduce the population A'^^-. (b) 
An ensemble of 7 x 10^ atoms are trapped within the TEMqo 
mode of an optical cavity using a ID intra-cavity optical lat- 
tice at 823 nm. Rotations of the Bloch vector are accom- 
plished using resonant microwaves at fa = 6.834 GHz. The 
photodiode records cavity transmission as the probe laser 
frequency is swept across the splitting, (c) The population 
= {Q-f/2g)^ is determined from the splitting Q,^ obtained 
by fitting the cavity transmission versus probe detuning to 
Lorentzians. (d) Experimental sequence. After preparing the 
initial Bloch vector Ji on the equator through optical pump- 
ing and a microwave 7r/2-pulse, a test rotation Rt is applied, 
and the result Jz — {Nf — Ni)/{N-f + Ni) is measured using 
the probe. 



is obtained from the measured and the J° for the same 
rotation without phase modulation using jz = Jz — J°- 
For future reference, jz corresponds to the special case 
of rotation axis about x, i.e. (pij — 0, and initial Bloch 
vector Ji = X. 

The nondemolition measurement is implemented by 
measuring the size of the collective vacuum Rabi splitting 
[41] generated by the coupling of a degenerate optical 
cavity mode to the |t) |e) transition, where the opti- 
cal excited state |e) = \F' = 1, mp = 0). The size of the 
splitting 51^ depends only on the total atomic population 
in It) as il^ = y/N^2g, where 2g is the single atom vac- 
uum Rabi frequency, a coupling constant determined by 
accurately known atomic properties and cavity geometry. 
The size of the splitting is measured by sweeping a probe 
laser across the resonances and fitting the transmitted 
power to Lorentzians. The population A'^ is determined 
from the measured splitting using A'^ = {fl^/2g)'^. Re- 
peating the same procedure after a microwave 7r-pulse 
swaps the atomic populations determines the population 
N^. From the measured populations A^^ and A^^, we ob- 



tain the quantity Jz 
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FIG. 3. Phase Modulated Rotations. Examples of the 
experimentally measured evolution of the Bloch vector along 
the measurement axis z versus rotation angle ip, or equiva- 
lently time, for (a) resonant modulation fm = fn = 40 A kHz, 
with /3 = 0.0125 rad, and (b) non-resonant modulation 
fm ~ 1.125/ij, with /3 = 0.025 rad. The modulation phase 
Om is for both (a) and (b). Dashed lines joins the data 
points (solid circles) to help guide the eye. The system re- 
sponds with frequency components at fm and leading to 
the observed amplitude modulation of the response in (b). In- 
sets show the ideal theoretical spiral trajectories of the Bloch 
vector about x. 



B. Response to Coherently Phase Modulated 
Rotation 

The harmonic oscillator-like response to a phase mod- 
ulated rotation is experimentally demonstrated in Fig. [3] 
The Bloch vector is prepared along x, then rotated nom- 
inally about X using microwaves whose phase is modu- 
lated at a fixed frequency fm- After a variable rotation 
angle the projection jz is measured. For phase modu- 
lation near resonance fm ~ fn , the envelope of the modu- 
lation grows roughly linearly with tp, whereas away from 
resonance fm 7^ fn, frequency components at fn and 
/,„ beat against one another to create amplitude mod- 
ulation. The insets of Fig. [3] show the ideal theoretical 



spiral trajectories about x for both cases. 



IV. NOISE IN BLOCH VECTOR DUE TO 
PHASE NOISE IN A SINGLE ROTATION 

A. Covariance Transfer Matrix 

Having obtained and experimentally demonstrated the 
response of the Bloch vector to a coherently phase mod- 
ulated rotation in Sec. |llj the goal of this section is to 
define a covariance transfer matrix that will allow the 
computation of the variance of the final Bloch vector 
projection along any arbitrary axis n due to a randomly 
phase modulated rotation caused by phase noise in the 
LO. 

We begin with modelling phase noise at a single dis- 
crete frequency /o by allowing the modulation phase 
and modulation amplitude /3 to take on random values 
between realizations of the rotation. Statistical results 
are obtained via ensemble averaging over all possible re- 
alizations of the phase modulated rotation. The modu- 
lation phase and amplitude are fixed in a single realiza- 
tion but random from one realization to the next. By 
drawing the modulation phase am from a uniform dis- 
tribution between and 2tt, and the modulation ampli- 
tude P from a Gaussian distribution with zero mean and 
variance we model the SSB phase noise C{fm) = 

{i3'^)6{fm - /o)/4 (see Appendix |b] for details). 

We define the covariance transfer matrix as the outer 
product T(0fl,J},V,/m) = 4(j/j})/(/?2), where (•) de- 
notes a statistical average over am and /3. The normal- 
ization of T (labels suppressed) is chosen so that inte- 
grating ■ T h over the SSB phase noise C{fm) of the 
LO yields the noise variance ((j/ • n)^) of the final Bloch 
vector projected along a i nea sure men t axis n. Using the 



solutions for jy,jz in Ec^. 10 and |ll[ we find the covari- 
ance transfer matrix T{ip, fm) for the special case i/);? = 
and J, = X 



Tyy{lp, fm) 

Tzz{lp, fm) 
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(l-a;2)2 
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(l-a;2)2 



(cos ip — COS (xtp))^ + {x sin tp — sin {xtjj))'^ 



(cos "0 — cos [xip))^ + (sin ip — x sin {xip)y 



(14) 

(15) 
(16) 
(17) 



Tyzi,"^, fm) = Tzyiip, fm) = z y (cos-0 - COS [x^)) sintp . 

1 — X'^ 

I 

The covariance transfer matrix T for arbitrary rotation tion = (J°, J°, J^) is derived using the small rotation 
axis azimuthal angle ipu and ideal final Bloch vector posi- r specified in Eq. [7] and [Sj The full analytic expression 
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FIG. 4. (color online). Transfer Functions. Examples 
of experimentally measured transfer function T^^i^l), fm) for 
tp — n (black filled circles), 2tt (red hollow circles) and 
47r (blue filled squares) overlaid on parameter-free theoret- 
ical curves. The amplitude of the phase modulation was 
P = 0.125, 0.05, 0.0625 rad for the i/- = tt, 2tv, i-jr transfer func- 
tions respectively. The modulation amplitude /3 ^ 1 is chosen 
to keep the response small to remain in the linear regime, and 
yet large enough to resolve the nulls at integer multiples of /jj. 
The transfer function does not depend on /3 as long as /3 ^ 1 
because both the response jz and the modulation amplitude 
are proportional to /3. 



for T is cumbersome but can be conveniently obtained 
through the transformation 

T = D{c^n, J})f{yj, fm)D(cj,R, J})^ , (18) 



where 



(Q -Jl COS c^R 

J° -Jlsin(j)R j . (19) 

^ J° cos -f J° sin (/>fl J 



As an example of this gene£al result, we experimen- 
tally measure the projection jz after rotating = x 
about X through different angles ip — tt, 27: and 47r 
while phase modulating the microwave source. Inte- 
ger multiples of tt were chosen to minimize sensitivity 
to intrinsic phase noise of the microwave source near 
DC. The transfer function T^z is obtained from aver- 
aging over four discrete values of the modulation phase 
oim — {0, f TT, ^}, while keeping the amplitude j3 con- 

stant, asfzzii>J,n) = ^ ELo i?(0, V', A /m, T)- The 
measured and theoretical transfer function Tzz are shown 
in Fig. |4| 

B. Covariance Noise Matrix 

In the linear response/small signal limit, the noise in 
the Bloch vector in some bandwidth is simply the integral 



of the noise variances due to phase noise at frequency fm 
over the relevant frequency bandwidth. Integrating the 
covariance transfer matrix T over the SSB phase noise 
^ifrn) of the LO yields the covariance noise matrix, de- 
fined for the largest possible bandwidth, 



F(^)= / T(^, /„)£(/„) d/„ 
Jo 



(20) 



The covariance noise matrix transforms to arbitrary (j)R 
and 3°j: in the same manner as the covariance transfer 
matrix in Eq. [18] via 

V = Dic^R,3})Vi^j)D{<j>R,J}y . (21) 

Using the covariance noise matrix, the variance in the 
projection along n may be obtained as 



((j/.n)2) ^nT.y.fi. 



(22) 



As a useful example of quantifying the noise mapping 
from the LO onto the Bloch vector, consider a white noise 
spectrum C{fm) — Co- Integrating the white phase noise 
spectrum over the covariance transfer matrix r(i/), fm) 
yields 

V{i!) = /loNEB , (23) 
where the noise equivalent bandwidth matrix is 



NEB = 7r/HSgn(V') 





i sin 2')/' 
■ sin^ ipi 





- sin^ 4) 

- i sin2V'y 



(24) 

A corollary to Eq. 24 is that as \tp\ increases, the covari- 
ance transfer matrix T becomes more and more sharply 
peaked around the Rabi frequency /r. Therefore for 
large 3> 1, most of the Bloch vector noise contribu- 
tion comes from phase noise near the Rabi frequency. As 
{ipl — >■ cxD, the covariance transfer matrix T approaches a 
delta function at /r 



Ti-tp, fm) ^ TrfR\^\Co S{fm - fR.) 




(25) 



The noise mapping for any rotation axis in the x — y 
plane can be obtained by transforming to arbitrary (/jr, 
but keeping = x. We find the non-zero elements are 

Vzz = Vzz COS'^ (j)R, Vyy = Vyy , Vy z = Vy z COS (j) R . 

Note that the variance Vzz can be driven to zero by ap- 
plying the rotation perpendicular to the Bloch vector as 
shown in Fig. [5jb). Alternately, keeping (pR = 0, but 
letting Jy = xcos6' + zsin0, i.e. the ideal final vectors 

lie in the x — z plane, one finds Vzz,yy = Vzz,yy cos^ 
and Vxx = Vzz sin^ 0. This noise mapping is graphically 
shown in Fig. [5ja). 
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FIG. 5. (color online). Noise Mapping, (a) Phase modula- 
tion of the LO causes modulation of the rotation axis (green 
arrow) about its mean orientation, here along x or (j>B, = 0. 
The final Bloch sphere of points is deflected by_an arnount de- 
scribed by a small rigid rotation r(0) = Ry{—jz)Rz{jy)- The 
deflection of a Bloch vector depends on the ideal final vector 
J^. If a set of possible deflections for an ideal final Bloch 
vector J°f along x is described by a circle jy + = const 
(shown outside sphere for clarity), the same set of deflections 
for other are described by ellipses and lines centered at 
J°f. (b) The noise mapping shown in (a) for (j)R — and for 
J J on the equator can be equivalently demonstrated by keep- 
ing the final vector oriented along = x and varying the 
rotation axis (pR. The observed noise variance Vzz (solid cir- 
cles) of the vector projection along z varies as the predicted 
cos'^ cj>R (solid line). The experiment was performed by apply- 
ing phase modulation at a discrete frequency fm ~ 1.125/_r 
to a rotation about x cos 4- y sin (pR that nominally rotates 
the Bloch vector = x cos 2<I)r + y sin 2(j}R through an angle 
^ = TT to = X. Averaging over the four modulation 
phases am ~ {0, f ,7r, ^} simulates phase noise at a single 
discrete frequency fm- 



V. NOISE IN BLOCH VECTOR DUE TO 
WHITE PHASE NOISE FROM MULTIPLE 
ROTATIONS 

A. Noise Propagation 

The single-rotation covariance noise matrix V allows 
us to analyze two crucial building blocks for coherent 
manipulation of quantum systems - composite pulses, 
used to suppress static amplitude and detuning errors, 
and pulse sequences designed to reduce qubit decoher- 
ence, such as spin echo and dynamical decoupling type 
sequences. Constituent rotations in composite pulses are 
applied in a back-to-back manner, leaving as little time 
as possible between the rotations. In contrast, time sepa- 
ration beween rotations in a pulse sequence may be com- 
parable or much longer than the time it takes execute a 
rotation. 

Assuming a white phase noise spectrum £(/m) = ^o, 
noise from different rotations become statistically in- 
dependent regardless of time separation between rota- 
tions. The time separation between rotations is, how- 
ever, important when dephasing of the Bloch vector sets 
a timescale Tc for the Bloch vector to phase diffuse by 



~ 1 rad. In that case, it becomes important to ensure 
that the shortest time separation r in the sequence satis- 
fies r ^ Tc so that the opening angle between the Bloch 
vector and rotation axis is well defined in order to ap- 
ply the formulae presented below. This is a reasonable 
condition given that pulse sequences designed to reduce 
decoherence also require r <C Tf, in order to be effective. 
In this context, the formalism for treating noise from a 
composite pulse, and noise from a pulse sequence are the 
same. 

We now present the noise propagation that gives the 
multiple-rotation covariance noise matrix W. The sym- 
bol W for the multiple-rotation covariance noise matrix 
is chosen to differentiate it from the single-rotation co- 
variance noise matrix V . A pulse sequence consists of N 
rotations with the fcth rotation given by 

i?fc = i?(0fe,Vfc). (26) 

In the absence of noise, the ideal Bloch vector after the 
fcth rotation is 

31 = Rk--- i?i J. . (27) 
The added noise from just the A:th rotation is 

Vk = DicPk,3l)V{i:k)D{^k,Jiy . (28) 

By accounting for how the noise from the previous ro- 
tation Wk-i is transformed by subsequent rotations, the 
total noise after the fcth rotation may be computed iter- 
atively using 

Wk = RkWk-iRl + Vk , (29) 

with the initial condition Wi = V\. Finally, the noise 
variance along any arbitrary projection axis fi after the 
fcth rotation is given by 

((jfc-n)2>=fiT.I^, .ft, (30) 
where jfc is the noise deflection after the fcth rotation. 

B. Average Infidelity 

While matrix elements of the covariance noise matrix 
W depend on pulse sequence specifics and the initial 
Bloch vector, we present here a simple formula that eval- 
uates the average quality of a pulse sequence using only 
general properties of the pulse sequence. Within the 
quantum control and computing community, the state 
infidelity 142] 

1 - F = Tr(VF)/4 (31) 

is an important measure of the rotation quality. While 
the state infidelity depends on pulse sequence details, the 
state infidelity averaged over the Bloch sphere of possible 
initial states 

(1 - F) = ^-TT/fl/Io/S (32) 
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depends only on the total rotation angle = X^fcLi V'fc 
and not the rotation axes (pk- Thus, a pulse sequence 
with smaller ^ is preferred over one with larger if 
one is mainly concerned with the average fidelity. In 
spin echo and dynamical decoupling schemes, suppression 
of environment-induced decoherence typically improves 
with the number of pulses. However, this comes at the 
expense of increasing the average phase-noise-induced de- 
coherence. It is, therefore, necessary to strike a balance 
between reducing environment-induced decoherence and 
reducing phase-noise-induced decoherence. 



C. Composite 7r-pulse Comparisons 

Composite pulses, designed to suppress static ampli- 
tude and detuning errors, have been thoroughly analyzed 
in the literature with regards to the degree of error can- 
cellation P^22j . The influence of phase noise on the 
Bloch vector through composite pulses, however, has re- 
ceived little attention in literature. 

Applying the noise propagation formalism presented 
in Sec. |V A| to the commonly used composite 7r-pulse se- 
quences: CORPSE, SCROFULOUS, and BBl JHHID], 
which effectively implement a 7r-pulse about x, we sum- 
marize in Fig. [6] the final variance Wzz and the state in- 
fidelity 1 — F versus the initial Bloch vector 3i specified 
by its polar angle 9i and azimuthal angle (f)i. Expres- 
sions for the composite pulse rotation sequences, covari- 
ance noise matrices, and state infidelities are provided in 
Appendix [Cj 

The single quadrature variance Wzz is of interest 
in metrology applications, particularly for manipulat- 
ing spin-squeezed states, as added noise in the squeezed 
quadrature can potentially destroy the squeezing. The 
state infidelity 1 — F, which includes variances from the 
two transverse spin components perpendicular to the 
ideal final Bloch vector, is particularly pertinent in quan- 
tum control for quantifying the overall quality of the ro- 
tations. 

It is generally not possible to minimize phase noise sen- 
sitivity and optimize static error cancellation simultane- 
ously. To understand the tradeoffs between phase noise 
sensitivity and static error cancellation, we compare and 
contrast the two for initial Bloch vectors in the x — y 
plane, i.e. Oi = 0, leaving (/j^ as the only degree of free- 
dom. We use Wzz and 1 — as the basis for evaluating 
sensitivity to phase noise. 

A static fractional amplitude error e results in an er- 
ror cTp in the rotation angle, and a static detuning error 
5 [fho- fa)/ Jr, where Jlo is the LO frequency, causes 
the rotation axis to be tilted up from the x — y plane by 
an angle arctan((5). The degree to which the static errors 
e, 5 are suppressed can be characterized by the static er- 
ror squared Wzz,st — j1 sti.^^ and the infidelity due 
to static error 1 - Fst ^ Tr(Wst(e, (5, (/)i))/4. These defi- 
nitions, in direct analogy to corresponding quantities for 
phase noise, allows meaningful comparison of the phase 
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FIG. 6. Composite 7r-pulses. Variance Wzz and infidelity 
1 - F of CORPSE TT-pulse (a, d), SCROFULOUS 7r-pulse (b, 
e) and BBl 7r-pulse (c, f) versus initial Bloch vector orienta- 
tion {9i,(j)i). 6i is measured from the x — y plane and 4>i is 
measured from x. Contours levels are normalized to JrCo- 
For scale, Wzz{di — 0,(j>i = 0) for a simple 7r-pulse about 
X is -K^ fuCo- Points of best (solid circle), worst (hollow cir- 
cle), and same (hatched circle) order of static amplitude error 
cancellation, and points of best (solid square), worst (hollow 
square), and same (hatched square) order of static detuning 
error cancellation (offset vertically for clarity) are shown for 
Bloch vectors with 6i — 0. All plots have the same axes as 
(a). A quadrant of the Bloch sphere is shown here. The rest 
of the contour plot can be generated using reflection symme- 
try about the x — y and x — z plane. The rest of the static 
error cancellation points can be generated via (f)i (j)i + tt. 



noise sensitivity and static error cancellation on the same 
footing in the next section. 

1. CORPSE -K -pulse 

The CORPSE 7r-pulse, used to suppress static de- 
tuning error, has the best static error cancellation at 

= Q,TT where both Wzz,st, 1 — ^st — 0{6^). However, 
it is also most sensitive to phase noise at (f)i = 0, tt (see 
Fig. ^sl, d)). At = 7r/2, the impact of phase noise is 
minimized at the expense of static error cancellation as 



both Wz 



1 - Fs, 



0{S^) + 0(s 



2. SCROFULOUS -K-pulse 

In contrast to the CORPSE 7r-pulse, sensitivity to 
phase noise does not vary with 0; for the SCROFULOUS 
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TT-pulse designed to suppress static amplitude error. As 
shown in Fig. [6]jb, e), one can simply choose (fii to op- 
timize the cancellation of static errors depending on the 
quantity Wzz,st or infidelity 1 — Fgt to be optimized, and 
on the dominant source of static errors (amplitude or de- 
tuning) without altering the impact from phase noise. 

3. BBl TV-pulse 

Finally, the BBl 7r-pulse, which compensates for am- 
plitude error with little to no cost to the sensitivity to de- 
tuning error, has similar impact from phase noise where 
the order of amplitude error cancellation is best (worst) 
at « 0.797r(0.297r) as shown in Fig. |6];c, f). There- 
fore one may choose to operate at « 0.797r where 
Wzz^st = 0{S^) + 0(ei") and 1 - i^st = 0{6^) + 0{e»). 
In fact, the impact from phase noise is slightly lower at 
4>i « 0.797r compared to at (j>i ~ 0.297r. On the other 
hand, there is a tradeoff between suppressing detuning 
errors and sensitivity to phase noise. Detuning error can- 
cellation is best at 0i = 7r/2 as both Wzz,st, 1 — ^st do not 
scale with S to any order if e = 0. However, the impact 
of phase noise is also worst at ipi = 7r/2. 

In general, careful evaluation of the relative scalings 
and contributions of phase noise, static amplitude and 
detuning errors is required to optimize the overall fidelity 
or specific noise quadratures. 



amplitude error as the Bloch vector nominally retraces its 
path while choice (a) accumulates static amplitude error 
as the Bloch vector keeps rotating about the same axis in 
the same sense. On the other hand, choice (a) does not 
accumulate static detuning error while choice (b) does. 
In future work, we will explore the possibility of engi- 
neering noise properties of spin echo pulse sequences. 

VI. CONCLUSIONS 

In summary, we have developed a general framework 
for analyzing the mapping of LO phase noise onto noise 
projections of a Bloch vector and extend the mapping 
to pulse sequences for the case of white LO phase noise. 
Detuning or transition frequency noise can be handled 
via the mapping /3 — A//™ where A is the frequency 
modulation amplitude. Results for special but important 
and illustrative cases are presented, which experimental- 
ists can readily utilize for estimation or design. Future 
work will extend the analysis to non-resonant excitation, 
more complex spin echo or dynamical decoupling pulse 
sequences, and include the effects of 1/f and higher or- 
der phase noise, where time separation between pulses 
may no longer be ignored, and noise correlations between 
pulses play an important role. 

VII. ACKNOWLEDGEMENTS 



D. Spin Echo Pulse Sequences 

Spin echo and dynamical decoupling sequences consti- 
tute another class of manipulation protocols in quantum 
control and computing, important for suppressing qubit 
decoherence, or, for instance, to undo probe-induced de- 
phasing IMTUl [13] . We analyze here spin echo sequences 
of the form [t - R{(j)i,Tr) - 2t - i?((/)2,7r) - 2t • • • - 
R{(j)N,n) — t] using the formalism developed in Sec. V A 



to find the covariance noise matrix W for the sequence. 

We consider the following two choices of rotation axes: 
(a) rotation axis always along x, or, (b) alternating be- 
tween X and —X with the first 7r-pulse applied along x. 
Note that the widely used Carr-Purcell [T] , Carr-Purcell- 
Meiboom-Gill [5], and Uhrig Dynamical Decoupling [i] 
sequences are special cases of (a) corresponding to spe- 
cific orientations of the Bloch vector with respect to the 
rotation axis x. 

Writing the initial Bloch vector as 3i = {J^., Jy, Jl), 
the covariance noise matrix for both choice (a) and (b) 
reads 



'1- j: 



Wn = Ntt^/rCo 



snJxJI J*^ 



S^j^jy 



(33) 



where sn = ( — l)'^^^-*. While the noise properties for 
the two choices are the same, their sensitivity to static er- 
rors are different. Choice (b) offers cancellation of static 
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Appendix A: Phase Noise Definitions 

The power spectral density of phase fluctuations 
S(j,{fm) of an oscillator is the mean squared phase fluctu- 
ations (A(/)(/m))^ at frequency offset from the carrier 
in a 1 Hz measurement bandwidth 



5^(/„) = A0(/™)VHz. 



(Al.l) 



The phase noise S^{fm) has units of rad^/Hz, and in- 
cludes contributions from both upper and lower noise 
sidebands at ±/m. 

One can measure S^{fm) by mixing the oscillator- 
under-test with a reference oscillator at the same fre- 
quency, and with much lower phase noise. The phase of 
the reference oscillator is chosen so that the mixer out- 
put v{t) is proportional to the relative phase difference 
(j){t) between the oscillator-under-test and the reference 
oscillator. The power spectral density S^{fm) is com- 
puted from the autocorrelation function of 4>{t) via the 
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Wiener-Khinchin theorem as 



{<l){t)<t){t + T))te-'^^f^^dT, (A1.2) 



where /,„ hes in the range (0,c»). 

The definition of the single sideband (SSB) phase noise 

C-{fm) is 



(A1.3) 



The units for C{Jm) are rad^/Hz. It is also commonly 
expressed in the form 10 logj^Q £(/,„), which has units of 
dB below the carrier in a 1 Hz bandwidth (dBc/Hz). The 
mean squared phase fluctuations observed in a Fourier 
frequency range from // to fh is given by 



(A1.4) 



For an oscillator whose amplitude noise is much lower 
than its phase noise, the SSB phase noise £(/m) is 
equivalent to the ratio of the power of a noise sideband 
PssB(/m) at frequency offset in a 1 Hz measurement 
bandwidth to the power in the carrier Pear as measured 
on a radio/microwave/optical frequency spectrum ana- 
lyzer 



P 



(A1.5) 



Appendix B: Relation between Mean Squared 
Modulation Amplitude and Single Sideband Phase 
Noise 



frequency /o, as being drawn from a random distribution 
of ara and /3. The modulation phase am is uniformly dis- 
tributed from to 27r, and the modulation amplitude /? is 
Gaussian distributed with zero mean and variance (/3^). 
The time-averaged squared phase modulation should be 
further averaged over the distribution for «,„ and /3 to 
yield the statistical phase fluctuations as 



2 



(A2.1) 



where we made use of the fact that averaging over time 
t has the same effect as averaging over phase a,„. Using 
Eq. [ALU [AO] and [A2J| we obtain the relation between 
the SSB phase noise C{fm) and mean squared modula- 
tion amplitude as 



-Sifm - /o) 



(A2.2) 



Appendix C: Composite and Single 7r-pulse: 
Covariance Noise Matrix and Infidelity 

The covariance noise matrix W and infidelity I — F for 
the CORPSE TT-pulse, SCROFULOUS 7r-pulse, BBl tt- 
pulse and a single 7r-pulse, which all effectively implement 
a TT-pulse about x, are given below. The results assume a 
white LO phase noise spectrum C{fm) = ^o- We specify 
the initial Bloch vector by its polar angle 9i, measured 
from the x — y plane, and its azimuthal angle 4>i measured 
from X so that Ji = (cos 0i cos ^i, cos 6'i sin 1/)^, sin 6'i). 



We establish here the connection between the mean 
squared phase modulation amplitude (/3^), used to nor- 
malize the covariance noise matrix T in Sec. [IV A[ and 
the SSB phase noise C{fm)- We model the phase mod- 
ulation (j){t) — /3sin(27r/ot -f a™), at fixed modulation 



1. CORPSE TT-pulse 

For the CORPSE TT-pulse sequence 
P(G, f)i?(TT, ^)P(0, ^) (time ordering right to 
left), used to suppress static detuning error, we have: 



1 - F 



^tt/t^/Io ((i37t - 3%/3) cos^ B., sin^ 0, 4- (^13tt + 3^3^ sin^ 

-nf^Co ^13tt — SVS^ cos^ 9i cos^ 0i 

^tt/j^/Io (i3tt + 3Vtj cos^ 0i cos^ 0, 

Wy, = ^ttJrLo (i3tt - 3V3) cos2 9, sin 20, 



— ^Tr/fl/Io (l3n + SVS^ sin26'i cos(/)i 
T^^fR^o (l37r - 3^cos26ii + (i3tt + 3^3) cos^ 6i ( 



(A3.1) 

(A3.2) 

(A3.3) 

(A3.4) 
(A3.5) 
(A3.6) 

(A3.7) 
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2. SCROFULOUS yr-pulse 



For the SCROFULOUS 7r-pulse sequence 
iT!(f ,7r) i?(^,7r) i? (f ,7r) (time ordering right to 
left), used to suppress static amphtude error, we have: 





3 

= -n'^fnCo (2 cos^ 6, sin^ 0^ + sin^ 9, 


) (A3.8) 


Wyy 


3 

= ■^Tr'^fR^o (2 cos^ 6i cos^ (f)i + sin^ 6 


) (A3.9) 


w,. 


3 

= i^TT^fnCo COS^ 0i 


(A3.10) 


W^y 


3 

= ^yx = cos^ 9i sm2(pi 


(A3.11) 


Wy. 


3 2 

= Wzy ^ -^'^ fn^o sin26'i sm^^ 


(A3.12) 




3 

= W:,x = -tt'^JrCo sin 20i cos 0^ 


(A3.13) 


1 - F 


lb 


(A3. 14) 



3. BBl TT-pulse 



The BBl TT-pulse sequence 

R{0,Tr) R{(j)fi, n) R{3(j)R, 2n) R{(j) a, tt) (time order- 
ing right to left), where = arccos(— |) w 104.5°, is 
used to compensate for amplitude error with little to no 
cost in the sensitivty to detuning error. The covariance 
noise matrix W and infidelity 1 — F for the BBl 7r-pulse 



are: 

W^^ = ^-it^IrCo (4cos2 sin2 4 + sin^ 6,) (A3.15) 

Wyy = ^tt^/rCo (4cos2 0i cos^ -I- 3sin2 0,) (A3.16) 

W,, = K^JrCo cos^ 0, (2 - cos 24>i) (A3.17) 

W^y = Wy^, = ^n^fRCo cos2 0, sin 2(j), (A3.18) 

Wy, = W,y = -^TT^fRLo sin 20, sine/., (A3.19) 
o 

W^, = W,^ = I^^JrC, sin 29, cos 0, (A3.20) 
o 

5 

1 - F = —n^fRCo (4 + 2 cos^ 61, - cos^ 9, cos 20,) 
16 

(A3.21) 

4. Single yr-pulse 

Finally, we provide expressions for the covariance noise 
matrix and infidelity for a single 7r-pulse R (0, tt) as a 
useful benchmark to compare against composite 7r-pulses. 
The noise covariance matrix W and infidelity 1 — _F for 
a single 7r-pulse are: 

W,, = tt^rCo (1 - cos2 0, cos" (A3.22) 
Wyy = TT^fRLo cos" 61, cos" (j), (A3.23) 
W^.. = ^'/iJ^o cos" 0, cos" (A3.24) 

W,y = W^y, = ^7r"/^/:o cos" 0, sin 20, (A3.25) 

^Wy,=0 (A3.26) 

VF,. = W'z. = ^7r"/fl£o sin 20, cos 0, (A3.27) 

1-F= ^TT^fRCo (1 + cos" 9, cos" (?!)i) (A3.28) 
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